Elimination of the Landau Ghost from Chiral 

Sohtons 

Josef Hartmann^ and Friedrich Beck 

Institut fuer Kernphysik 
Technical University of Darmstadt 
SchloBgartenstr. 9, D-64289 Darmstadt, Germany 

and 

Wolfgang Bentz 
Department of Physics, 
Faculty of Science 
University of Tokyo 
Kongo 7-3-1, Bunkyo-ku, Tokyo 113, Japan 

Abstract 

We show a practical way based on the Kanen-Lehmann represen- 
tation for the two-point functions to eUminate the instabihty of the 
vacuum against formation of smah sized meson configurations in the 
chiral a model. 



Supported by Deutsche Forschungsgemeinschaft under contract BE 348/11-1,2 



1 Introduction 



Since QCD is presently intractable at low energies, approximation schemes 
based on effective models have been introduced |1|]. These models in a way 
incorporate some important features of QCD, like chiral symmetry and its 
spontaneous breakdown. 

Though being approximative, these effective models still cannot be treated 
exactly. Therefore we have essentially two levels of approximation: i) when 
we choose the model itself and ii) when we specify an approximation scheme 
in order to perform actual calculations within the model. This implies that 
the results may depend on the actual approximation scheme. Therefore, a 
failure of the model could be due to the model itself, due to the approx- 
imation scheme or both. So, if we want to test the model, we should be 
careful about the approximation scheme and keep the dependence on ad hoc 
parameters as small as possible. 

One of the effective models of QCD to describe the nucleon is the chiral 
(T-model [§. The standard procedure to calculate nucleonic properties in the 
soliton picture is to use the loop expansion up to the one quark loop level, i.e; 
to integrate out the quark degrees of freedom and use the remaining purely 
mesonic effective action at the zero loop level in the mesonic fields 0. One 
can choose some cut-off procedure to regularize the quark determinant [^], 
but in any case one has to renormalize by introducing appropriate counter 
terms in order to define the parameters appearing in the original Lagrangian 
of the model p. Eventually one can let the cut-off go to infinity as it is 
usually done in renormalization theory, or leave it finite. We note that in the 
(T-model, however, neither the cut-off scheme itself nor the size of the cut-off 
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can be determined uniquely. 

It is known, however, that the loop expansion for asymptotically non free 
theories without ad hoc cut-off parameters displays an unphysical pole in 
the corresponding propagators, the so called Landau ghost Therefore, 
if one follows the procedure outlined above, one finds that the renormalized 
(T- model is plagued by the Landau ghost, which leads to an instability of the 
usual translational invariant vacuum 0, In view of our remarks at the 
beginning of this section, the problem of this Landau ghost can be treated 
in several ways: 

1. The effective model is considered as unphysical and is abandoned. 

2. The approximation scheme (loop expansion) is considered appropriate, 
but the effective model can be used only with a regularization scheme, 
hereby being restricted to energies far below the unphysical pole. 

3. The approximation scheme is improved such that the effective model 
does not a priori yield unphysical results for energies up to the order 
of the baryon mass. 

The second viewpoint implies that vacuum loops are either simply discarded, 
or calculated with finite cut-off. Concerning the first alternative (drop vac- 
uum loops), we note that in some situations, like in an external field, such a 
treatment violates conservation laws and cannot be implemented consis- 
tently. Besides this, the investigation of the change of the vacuum structure 
due to the finite density should be an important subject for any relativistic 
field theory. Concerning the second alternative (finite cut-off), we note that 
the Landau ghost in hardonic theories occurs at rather low energies of the 
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order of the baryon mass (~ IGeV). This means that its presence is a real 
problem for effective quark theories, and its avoidance by the introduction 
of an ad hoc parameter is somewhat unsatisfactory. Since the energy scale 
relevant for effective quark theories is set by hadronic masses, a treatment 
according to the third viewpoint, which we will adopt in this paper, is a fea- 
sible alternative. We consider the Landau ghost as one of the basic infinities 
in effective field theories which should be brought under control without the 
introduction of further parameters. In a different context, this viewpoint 
has been taken already a long time ago by Redmond and Bogoliubov et al 
[p!0[] who have shown how to construct ghost free propagators in the frame- 
work of the loop expansion based on the requirement of the Kallen-Lehmann 
representation. Their propagators have several interesting features: The as- 
sociated wave function renormalization constant is finite, and they have an 
essential singularity at = {g is the coupling constant) which is expected 
from more intuitive physical arguments [0. Moreover, since the associated 
running coupling constant has no singularity at space-like momentum trans- 
fers, the 'zero charge problem' is avoided. Their method has recently been 
formulated in the language of the effective action |T^, which is suitable for 
effective quark theories, and also applied to infinite systems [|^, [l^. The 
main purpose of the present paper is to show a practical way to incorporate 
the ghost elimination in a finite system like the soliton, and to demonstrate 
that with this procedure the translational invariant vacuum becomes stable 
with respect to decay into an array of small sized configurations. We will not 
construct the fully self consistent solutions including valence quarks in this 
paper which is left for a future work. (For another possibility to avoid the 
Landau ghost by introducing vector bosons, see ref. IQ.) 
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The cr-model has so far not yet been fully exploited to obtain self consis- 
tent solutions including the effect of the Dirac sea. One of our motivations 
for investigating this model is as follows: In contrast to the Nambu-Jona- 
Lasinio (NJL) model ||T^, which has been extensively used to construct soli- 
tonic solutions |jl6|, |1^, |18[ , it has a 0^ interaction term with variable strength 
(specified by the a mass) which controls the deviation of the chiral radius 
from its vacuum value. Such a term would contribute to the ground state 
energy the more the farther the chiral fields are away from the chiral circle 



of the vacuum. Since it is known that interactions terms of the form 6^ [119 



or of higher order like the t'Hooft determinant interaction [T^ can prevent 
the collapse of the NJL soliton, there is the possibility that stable solitons 
exist for the linear cx-model. The major obstacle to a full investigation of 
the solitonic sector of the cx-model, however, was the instability of the trans- 
lational invariant vacuum. The fact, to be established in this paper, that 
the method of Redmond and Bogoliubov et al leads to a stabilization of the 
translational invariant vacuum is a first step towards self consistent solutions 
in the cr-model. 

The rest of the paper is organized as follows: To be self contained and to 
set the formalism we state in sect. 2 some results for the ground state energy 
in the cr-model obtained earlier |^ , and explain the role of the Landau ghost 
in a way which seems most transparent to us and which is appropriate for 
implementing the ghost elimination procedure. In sect. 3 we construct a 
ghost free model and show results for the ground state energy, and in sect. 
4 we summarize our results. 
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2 The ground state energy and the Landau 
ghost 



We now proceed to exhibit how the Landau ghost contributes to the vacuum 
instabihty. The Euchdean Lagrangian of the cr-model reads, after introducing 
the 'physical' parameters rricr and m,r instead of the original ones /i^ and 



with 



M 



C = Cf + Cm + Csb (1) 

Cf = ^Dtjj; D = + MU (2) 

id,U) {d,U+) + \{ml- ml) {U+U - if + ml (U+U - 1)1 . (3) 



Here ip stands for the u and d quarks, M = gv is the effective quark mass, 
U = -(cr + ^7577 • r) is the chiral field, and Csb = fnfnl cr with f.^ the pion 
decay constant. The model contains the two free parameters g (or M) and 
mcr- V is the vacuum expectation value of cr, and in the symmetric limit 
(m^r = 0) we have v = fn^- Our Euclidean metric is such that a;^ = 
= (r, r) with r = < t < j3 with j3 the upper limit of Euclidean 
time integration, and 7^ = 7^ = {iPD,j) with the usual Dirac /^-matrix. 
Defining the effective bosonic action F from the generating functional Z as 
usual by Z = / DU exp(— F) we obtain after integrating out the quark 
fields 

T = Tf + Tm + Tsb (4) 

with 

Tp = -\Tr In D+D + c.t. = -\Tr ln(l + GV) + \TrGV - {TrC^YXb) 
Tm = J (I'^xCm; TsB = Id^xCsB- (6) 
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Here TtA = Nc J d'^xtT{x\A\x) where Nc = 3 is the number of colors and tr 
refers to the Dirac and isospin indices, G = (—9^ + M^)~^, and 

V = iM^U + M^{U+U -1). (7) 

The integrals over the Euclidean time xq extend from to f3. The counter 
terms (c.t.) in which include the subtraction of the translational in- 
variant vacuum contribution {U = 1), are determined such that Tp gives no 
terms of the same form as those already present in Tm- This corresponds 
to the mesonic mass and wave function renormalization at the renormaliza- 
tion point /i^ = f^. From here on we consider the chiral symmetric case 
(m^ = 0). 

The vacuum instability can be seen by expanding the ground state energy 
E = V / (3 = (Tp + Tm)/i3 in powers of a characteristic length scale R where 
the classical meson fields are localized. We assume time independent fields 
f/(r) = ?7(x) with x = r/i?. If we use x as the integration variable, every 
derivative gives rise to a factor 1/R, and therefore terms involving deriva- 
tives will be dominant for small R compared to those without derivatives. 
Therefore, up to 0{R) Tm contributes only the kinetic term. To obtain the 
lowest order contribution from Tp we expand the logarithm in (^) and note 
that the term linear in V is cancelled by a counter term. Then for small R 
the leading term is the one quadratic in the derivatives, i.e; quadratic in the 
first term of eq. (0). From the above we see that the linear and the non 
linear model look the same at small sizes. The complete term of order is 

T^;-'^ = hrr^GVGV - (8) 

Let us consider eq. (H) in a continuous and infinite plane wave (PW) basis 
{x\k) = exp(i/cx), where is the volume. Then the matrix element 
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(A;'|V^|A;) equals Q^^ 6qQflV{q), i.e; the Fourier transform of V{r) with q = 
k — k'. Using the dimensionless variables x = r/i? and t = qi? and keeping 
for V only the first term in eq.(|^) we obtain for (P) in the continuum limit 



lTr^^\NcM'Rj0^t.\mu[QyO(R') (9) 

with 

'^^^'^ " / (2^ ( (A;2 + M2) ((A; + q)^ + M^) ~ + M2)2 ) 



and L(t) is the Fourier transform of (^^^t/(x)j: 

L(t) = j e** " ^^U{x)£x. (11) 

Adding the kinetic term, we get the leading contribution to the effective 
action in the small size expansion as 

1 q^ 

Since — ^ In — - as — > oo it follows that (Bl) behaves as 

167r2 M2 "-^ 

a M^R \n{MR) for i? ^ with a > 0, which overcomes the kinetic term (the 
1 in ((|1^))) for R small enough, leading to a negative value for r*^^ ''\ This 
means that for small sizes R the energy of this 'non-translational invariant 
vacuum' can become lower than the energy of the translational invariant 
vacuum. 

To see that the Landau ghost in the meson propagators is responsible for 
this vacuum instability, we use the fact that the effective action, when ex- 
panded around the fields in the translational invariant vacuum (a = f , tt = 
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0), gives the inverse a and vr propagators in the translational invariant vac- 
uum as the coefficients of the terms of second order in s = a — v and tt. If 
we express V of eq. (0) in terms of cr, tt using U = -(a + ^7577 ■ r), we see 
that these second order terms are completely contained in (^. Taking also 



the terms of second order in s, tt in the mesonic term in (pf), we obtain 13 



l^^'^ = \l ^M'^y<-^)G~^'^'^") + ^(q) • ^(-q)G';^(q^)) (13) 

with the inverse Euclidean Schwinger-Dyson propagators 

G-\q') = (1 + ViVc0(g2)) , (14) 
G-\q') = q'{l + 4g^Nc<P{q'))+ml + 16g^NcM^(f){q'). (15) 

In fig. |l] we show the inverse pion propagator (|T1|) for g = 4 hj the dashed 
line. The previously noted behaviour of (plq"^) leads to the Landau ghost 
pole at Euclidean g^, and as a consequence (|T3|) can become more and more 
negative for large q^, i.e; for small sized fields. In this limit of small sized 
fields the leading contribution to (|13]) comes from the terms oc q^ in eqs. 



(|n|) and (|T3|), and if expressed in terms of t = qi? this gives again eq. (|^). 

8 1 
[Note that tr|L(t)P = — (s(t)5(-t) + n(t) ■ n(-t)) with §(t) = ^s(q) 

the Fourier transform of s(x) and similar for Trft).! 



3 Ghost removal in a finite system 

We now address the problem of removing the Landau ghost in a finite system. 
The procedure is described in ref. |1^ and has been extended to the path 



integral formalism in ref. [|T^. To calculate the ground state energy we 
have to ffist obtain meson propagators with the correct analytical properties. 
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These propagators are then implanted into the effective action which is in 
turn used to calculate the ground state energy. 



The prescription of ref. [T^ to eliminate the Landau ghost from the effec- 
tive action is to replace the Schwinger-Dyson propagators in the second 
order term (0) by the Kallen-Lehmann (KL) propagators Aq,. According to 



ref. |]T0[, these are constructed from the KL representation using a spectral 
function obtained from the one-loop meson self energy, and are free of the 
Landau ghost. Thus, in this method the loop approximation is is used only 
for the imaginary part, while the real part is calculated from the dispersion 
relation. Since this method avoids the Landau ghost, it clearly represents 
an improvement of the straight forward loop approximation for the whole 
propagator. The new effective action becomes 

Tkl = t + 5v, 5r = rg)^-r(2). (le) 

In order to be consistent with chiral symmetry, however, this ghost subtrac- 
tion has to be done under the constraint of the Ward identity p[ ^^{(f)^^ — 
A^(g^)^^ = ivT{~q; q,0) with T(— 0) the crvr^ vertex where one of the 
external pions has zero momentum. To preserve this identity the difference 
of the a and vr inverse propagators has to remain invariant under ghost sub- 
traction [|l3l, i.e; A^^ — = A^^ — G~^, and we obtain from (|16]) and 



^^r = i / ^(^(q)^(-q) + ^(q) ■ ^(-q))(A;^(q^) - G-\cf)). (17) 

It can be shown ||10| that the KL propagator can be obtained from the SD 
propagator by subtracting the ghost pole: A.j.^q'^) = G-^i^q^) — Zojic^ — m|.), 
where < is the residue of the pole at Euclidean = m^. We show A^^ 
for m,r = and = 4 by the solid line in fig. |1|. 
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The ghost free ground state energy is thus obtained by adding the piece 
(\n\} to -^r with r given by eq. (^). The exact evaluation of Tp, given by 



eq. (P), requires the diagonahzation of the Dirac hamiltonian h = —icx ■ V + 
PdMU, and this is done most conveniently in a basis |A) which diagonalizes 
the free Dirac hamiltonian ho = h{U = 1) = —ia ■ V + (3dM\ ho\\) = e°|A) 



with el = ykl + in a box of size D with discrete momenta Ik^l < 
kmax- (Here A labels all necessary quantum numbers except colour.) D and 
kmax should bc taken large enough such that all results are unchanged under 
further increase. We will call this discrete and finite basis the Kahana-Ripka 
(KR) basis ||20|. On the other hand, the ghost subtraction is formulated most 
conveniently in the continuous and infinite PW basis, for which the result 
is given by eq. (0). Therefore, we proceed in two steps: First, in order 
to demonstrate that the KR basis gives the same results as the PW basis 
provided that k^ax and D are taken large enough, we evaluate Fi? in the 
small size approximation both in the KR and the PW basis and compare 
the results. Then we evaluate Tp exactly in the KR basis and add the ghost 
subtraction term (p!7|). 

The expression for F^? in the small size approximation using the PW basis 
has already been given in eq. @, and the corresponding result for the KR 
basis is obtained from eq. (H) as 

Since for the exact calculation one has to diagonalize h, it is convenient 
here to treat (A'|y|A) as h"^ — where /i^^ = ha-yh^/s is the square of a finite 
dimensional matrix. In fig. |^ we show the small-size approximation to the 
Fermion loop energy, based on eq. (y) (solid line) and eq. (|18]) (dashed line) 
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in units of MNc as a function of Mi?, assuming a Hedgehog profile with the 
vacuum value for the chiral radius and winding number n = 1: 

U = expiii ■ TQ{r)'~^^), Q{r) = iTexp{—r / R). (19) 

(In this case there are no additional 0{R'^) terms in (P) due to U^U = 1. 
Also, Tf/M or r^^-'V^, viewed function of MR, is independent of M 
(or g). When the KR basis is used, the effective action becomes a function 
of kmax, and in order to reach convergence {T{kmax) T) in the small-size 
approximation we had to go up as high as kmax — 40M for MR ^ 1 . This is in 
contrast to models with finite cut-off like the NJL model, where kmax — lOM 
is sufficient for solitons of normal size ||T^. The calculations in this paper 
are performed with values for {kmax/M, DM) ranging from (70,5.7) for small 
R to (20,20) for large R, such that kmaxD ~ 400 is kept constant. Fig. || 
demonstrates that the results using the PW and the KR basis agree well if 
kmax and also the box size D are taken large enough. 

Comparison with the full calculation, to be discussed below, shows that 
the result of fig. | is quite accurate up to MR ~ 0.7. The fact that the 
ground state energy relative to the translational invariant vacuum is negative 
for these small R indicates the vacuum instability due to the Landau ghost 
discussed above. 

The total ground state energy in the small size approximation is shown in 
fig. ^. It is obtained from eq. (0), and for the ghost free model by further 
adding the piece (|17[). We see that the ground state energy in the small 
size approximation becomes a positive quantity after removal of the Landau 
ghost. 

We now return to the full effective action, eqs. (|§),(|§). In the KR basis 
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Tp can be written in the form 

with V = h"^ — a.s before, and e\ are the eigenvalues of h. Since we have 
already demonstrated for the small size expansion that for kmax and D large 
enough the results with the KR basis agree well with those of the PW basis, 
it is feasible to use the expression in the PW basis, eq. (]T7|), for the ghost 
subtraction term ST. 

Fig. |] shows the total ground state energy -^F including the ghost (dashed 
line) and ^^kl without the ghost (solid line) in units of MNc as a function of 
RM for the Hedgehog meson profile ([T9|). After ghost subtraction, the total 
ground state energy relative to the translational invariant vacuum is positive, 
which demonstrates that amending the two-point functions such that they 
satisfy the KL representation is sufficient to stabilize the vacuum against the 
formation of arrays of strongly localized meson field configurations. 

4 Summary 

In effective quark-meson theories, which have been designed to model QCD 
in the energy region of hadron masses, the Landau ghost appears in the me- 
son propagators at rather low energies of ~ IGeV. Its presence is a real 
problem for these theories and leads to an instability of the translational in- 
variant vacuum. In this paper we have taken the viewpoint that the Landau 
ghost should be brought under control without the ad hoc introduction of 
further parameters. Based on the work of Redmond and Bogoliubov et al 
1^ , who have shown how one can improve the loop expansion such that the 
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propagators have the correct analytical properties, we have shown a practical 
way to eliminate the Landau ghost from chiral solitons in the cr-model. We 



have demonstrated that the ghost free effective action formalism of ref. [12 



can be applied successfully to finite systems. In particular, since the ghost 
subtraction term is expressed most conveniently in a continuous and infinite 
plane wave basis, we have investigated in detail the conditions under which 
it is feasible to add this term to the Fermion loop term which is evaluated in 
a discrete and finite basis. By considering vacuum configurations character- 
ized by various sizes of the meson profiles we have shown that in our ghost 
free model the translational invariant vacuum has the lowest energy and is 
therefore stable with respect to decay into small sized configurations. 

The method given in this paper can form the base for further investiga- 
tions on self consistent soliton solutions in the cr-model. In particular, as 
mentioned in the Introduction, it would be interesting to see whether self 
consistent solutions can be obtained with a reasonable strength of the (j)'^ 
term. 

The problem of the Landau ghost and the associated vacuum instability 
occurs also in relativistic meson-nucleon theories for nuclear structure [PT|] . 
For example, if a finite nucleus is described in the relativistic Hartree ap- 
proximation including the effect of the Dirac sea, one has to apply the same 
'overall ghost subtraction' as performed in this paper (see eq. (|TB|)). In more 
sophisticated approximations like the '1/N' expansion p2| or the 'modified 
loop expansion' PSI the Landau ghost appears also in the subgraphs (even 



for infinite systems), and the ghost subtraction must be applied to these sub- 
graphs. Thus, although we used the language of effective quark theories in 
this paper, the method is applicable to relativistic nuclear structure physics 
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as well. 
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Figure caption 

1. The inverse pion propagator with m^r = 0, (7 = 4 for Euchdean q^. The 
dashed hne shows the Schwinger-Dyson propagator and the sohd hne 
the Kallen-Lehmann propagator. 

2. The Fermion loop energy Ef = ITf in units of MNc function of 
MR in the small size approximation. The solid line shows the result 
with the PW basis, eq.(^), and the dashed line shows the result with 
the KR basis, eq. (|TB|) with values for k^ax and D described in the 
main text. 

3. The total energy E = m units of M function of MR. The 
dotted and dashed-dotted lines refer to the small size expansion be- 
fore and after the ghost subtraction, respectively, and the dashed and 
solid lines show the full result before and after the ghost subtraction, 
respectively. 
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